CONVERGENCE OF A MASS CONSERVING ALLEN-CAHN EQUATION 
WHOSE LAGRANGE MULTIPLIER IS NONLOCAL AND LOCAL 



MATTHIEU ALFARO AND PIERRE ALIFRANGIS 



ABSTRACT. We consider the mass conserving Allen-Cahn equation proposed in (8): the 
Lagrange multiplier which ensures the conservation of the mass contains not only nonlo- 
cal but also local effects (in contrast with 1141 ). As a parameter related to the thickness 
of a diffuse internal layer tends to zero, we perform formal asymptotic expansions of the 
solutions. Then, equipped with these approximate solutions, we rigorously prove the con- 
vergence to the volume preserving mean curvature flow, under the assumption that classical 
solutions of the latter exist. This requires a precise analysis of the error between the actual 
and the approximate Lagrange multipliers. 



1. Introduction 

Setting of the problem. In this paper, we consider u e = u e (x,t) the solutions of an 
Allen-Cahn equation with conservation of the mass proposed in |8 |, namely 



(1.1) 3 t u e = Au e + — ( f(ue) - - Jn ; jT „7 = V^Wju,) ) inQx (0,oo), 



supplemented with the homogeneous Neumann boundary conditions 

on 

(1.2) —^(x,t)=0 ondft x (0,oo), 
and the initial conditions 

(1.3) u e (x, 0) = g e (x) in £1 

Here 51 is a smooth bounded domain in M. N (N > 2) and v is the Euclidian unit normal 
vector exterior to dil. The small parameter e > is related to the thickness of a diffuse 
interfacial layer. The term 

j n f(u e (x,t))dx 



(1.4) Jn , ex = = — JAW(u £ (x,t)) 

J n ^W(u £ (x,t))dx V 

can be understood as a Lagrange multiplier for the mass constraint 

d f 

(1.5) — / u 6 {x 7 t)dx = 0. 

ut Jq 

Let us notice that ( 11.4b combines nonlocal and local effects (see below). 

The nonlinearity is given by f(u) := —W'(u), where W(u) is a double-well potential 
with equal well-depth, taking its global minimum value at u = ±1. More precisely we 
assume that / is C 2 and has exactly three zeros — 1 < < +1 such that 

(1.6) /'(±1)<0, /'(0)>0 (bistable nonlinearity), 
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and that / is odd, which implies in turn that 

(1.7) / f(u)du = Q (balanced nonlinearity). 

The condition ( 11.61 ) implies that the potential W(u) attains its local minima at u = ±1, and 
( 11.7b implies that W{— 1) = M /r (+l). In other words, the two stable zeros of /, namely 
±1, have "balanced" stability. For the sake of clarity, in the sequel we restrict ourselves to 
the case where 

(1.8) f(u) =u(l-u 2 ), W(u) = \(l-u 2 ) 2 - 

This will slightly simplify the presentation of the asymptotic expansions in Section [4] and 
is enough to capture all the features of the problem. 

The initial data g e are well-prepared in the sense that they already have sharp transition 
layers whose profile depends on e. The precise assumptions on g e will appear in Theorem 
For the moment, it is enough to note that 



. m „ — 1 in the region enclosed by r 

(1.9) lim g s = < 

e^o 1+1 in the region enclosed between and To, 

where Tq CC SI is a given smooth bounded hypersurface without boundary. 

Our goal is to investigate the behavior of the solutions u e of ( 11.11 ). ( 11.21 ). ( 11.31 ), as e — >• 0. 

Related works and comments. It is long known that, even for not well-prepared initial 
data, the sharp interface limit of the Allen-Cahn equation dtu £ = Au e + e~ 2 f(u e ) moves 
by its mean curvature. As long as the classical motion by mean curvature exists, it was 
proved in lfT2l and an optimal estimate of the thickness of the transition layers was provided 
in . Let us also mention that, recently, the first term of the actual profile of the layers 
was identified Q. If the mean curvature flow develops singularities in finite time, then a 
generalized motion can be defined via level-set methods and viscosity solutions, Ifl8l and 
|[T5l . In this framework, the convergence of the Allen-Cahn equation to generalized motion 
by mean curvature was proved by Evans, Soner and Souganidis ifPTl and a convergence rate 
was obtained in Q). 

The above results rely on the construction of efficient sub- and super-solutions. Nev- 
ertheless, when comparison principle does not hold, a different method exists for well- 
prepared initial data. It was used e.g. by Mottoni and Schatzman l24l for the Allen-Cahn 
equation (without using the comparison principle!); Alikakos, Bates and Chen Q for the 
convergence of the Cahn-Hilliard equation 

(1.10) d t u e + A ^A« E + i/(« e )J = 0, 

to the Hele-Shaw problem; Caginalp and Chen ifTUl for the phase field system... The idea 
is to first construct solutions u e ^ of an approximate problem thanks to matched asymptotic 
expansions. Next, using the lower bound of a linearized operator around such constructed 
solutions, an estimate of the error ||u s — u e ^\\lp is obtained for some p > 2. 

Using these technics, Chen, Hilhorst and Logak 1141 considered the Allen-Cahn equa- 
tion with conservation of the mass 

(l.H) d t u e = Au E + i (f{u s ) -W\J Q f^H ' 

proposed by ll25l as a model for phase separation in binary mixture. They proved its 
convergence to the volume preserving mean curvature flow 

(1.12) V n = -K+^ 1 [ ndH"- 1 onr f . 

I 1 1\ Jr t 
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Here V n denotes the velocity of each point of T t in the normal exterior direction and k the 
sum of the principal curvatures, i.e. N — 1 times the mean curvature. For related results, 
we also refer the reader to the works [0 (radial case, energy estimates) and ll22ll (case of a 
system). 

In a recent work, Brassel and Bretin [8] proposed the mass conserving Allen-Cahn 
equation ( II. lb as an approximation for mean curvature flow with conservation of the vol- 
ume ( 11.12b . According to their formal approach and numerical computations, it seems that 
"( 11.11 ) has better volume preservation properties than (11.1 11 1". Let us notice that, as far as 
the local Allen-Cahn equation is concerned, such an improvement of the accuracy of phase 
field solutions, thanks to an adequate perturbation term, was already performed in 11201 or 
infl31. 

In the present paper we prove the convergence of (11.1) to (11.12b . Observe that in (11.11b 
the conservation of the mass ( 11.5b is ensured by the Lagrange multiplier — tjjt J*q f(u e ) 
which is nonlocal, whereas in the considered equation (11.1b the Lagrange multiplier (11.4b 
combines nonlocal and local effects. On the one hand, this will make the outer expansions 
completely independent of the inner ones, and will cancel the e order terms of all expan- 
sions (see Section|Ui. On the other hand, this makes the proof of Theorem l2.3l much more 
delicate since further accurate estimates are needed (see subsection l6.1b . In other words, 
in the study [ 14 1 of ( 11.11b . it turns out that the nonlocal Lagrange multipliers "disappear" 
while estimating the error estimate u e — u e ^- This will not happen in our context and our 
key point will be the following. Roughly speaking, our estimates of subsection 16. II will 
make appear an integral of the error on the limit hypersurface which must be compared 
with the L 2 norm of the error. If the former is small compared with the latter then the 
Gronwall's lemma is enough. If, as expected, the error concentrates so that the former 
becomes large compared with the latter, then the situation is favorable: a "sign minus" 
intends at decreasing the L 2 norm of the error (see subsection 16. II and Remark l6T2l for 
details). 

To conclude let us mention the work of Golovaty lED . where a related equation with a 
nonlocal/local Lagrange multiplier is considered. The convergence to a weak (via viscosity 
solutions) volume preserving motion by mean curvature is proved via energy estimates. 
As mentioned before, our method is different and allows to capture a fine error estimate 
between the actual solutions and the constructed approximate solutions. 



2. Statement of the results 

The flow ( 11.121 ). Let us first recall a few interesting features of the averaged mean curva- 
ture flow (11.12b . It is volume preserving, area shrinking and every Euclidian sphere is en 
equilibrium. The local in time well posedness in a classical framework is well understood 
(see Lemma |2~T1 for a statement which is sufficient for our purpose). It is also known that 
local classical solutions with convex initial data turn out to be global. Additionally, there 
exist non-convex hypersurfaces (close to spheres) whose flow is global. For more details 
on the averaged mean curvature flow (11.12b . we refer the reader to |[T9l . [23), lTT6l and the 
references therein. 

Lemma 2.1 (Volume preserving mean curvature flow). Let Slo CC Vibe a subdomain such 
that To := ddo is a smooth hypersurface without boundary. Then there is T max € (0, oo] 
such that the averaged mean curvature flow ( 11.12b , starting from Tq, has a unique smooth 
solution U <t<T m »-(r t x {t}) such thatT t CC Q, for all t <G [0,T max ). 

In the sequel, for T as in (O, we fix < T < T max and work on [0, T]. We define 



r := u < t <T(r t x {*}), 
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and denote by O t the region enclosed by T t . Let us define the step function u = u(x, t) by 

(2.1) u{x,t) := \ ~\ for all t e [0,T], 

I +1 in 12 \ S2t 

which represents the sharp interface limit of u e as e — > 0. Let d be the signed distance 
function to T defined by 



(2.2) d(x, t) 



-dist(x,r t ) for a; € Vt t 
dist(a;,r t ) foxxe£l\Q. t . 



Main results. We rewrite equation (ll.lt as 

(2.3) d t u £ ~Au £ -^ (/(« e ) - £ A £ (i) V4VK(u £ )) = inOx(0,oo), 
by defining 

(2.4) £ A £ (t) 



u F — U F 



Our first main result consists in constructing accurate approximate solutions. 

Theorem 2.2 (Approximate solutions). Let us fix an arbitrary integer k > max(iV, 4). 
Then there exist (u £tk {x, t), X £j k{t))xeQ o<t<T sucn that 



^ (f(u e ,k) ~ e\ e , k (t)^4W{u e , k ^j = S £ . k in fi x (0, T), 



(2.5) d t u £ ^ k - Au Ejfc 
with 

(2.6) UMU~(nx(0,T)) =0{e k ) as e -> 0, 
ana? 

(2.7) _^!*( x t) = o ondrix(0,T), 

ov 

(2.8) — / u £ . k (x,t)dx = for all te (0,T). 
»fc J si 

Observe that by integrating (12.5) over 17 and using (12.71 i and (12.8b . we see that 

/n V4^(u e ,fc) 

Then we prove the following estimate, in the L 2 norm, on the error between the approx- 
imate solutions u e> k and the solutions u £ . 

Theorem 2.3 (Error estimate). Let us fix an arbitrary integer k > max( TV, 4). Let u £ be 
the solution of (11.1b , ( 11.2b , ( ll.3l l wirfi f/ie initial conditions satisfying (II. 9b ant/ 

(2.10) ffe (x)=tl 8lfc (s,O) + e (s)G[-l,l], /<^e = 0, ||0 £ || L 2 ( n) =0(£ fe -^). 

77ie«, f/iere w C > rac/i f/zaf, /or £ > small enough, 



sup |K(-,t) - u £ , k (-,t)\\ L 2 iQ) < Ce k 

0<t<T 

Let us notice that, since —1 < g £ < 1, it follows from the maximum principle that 
— 1 < u e < 1. Also since g £ ^ 1 and g e ^ — 1, the conservation of the mass implies 
w e ^ 1 and u e ^ —1, which shows that the definition of e\ £ (t) in (12.4b actually makes 
sense. 

As it will be clear from our construction in Section the approximate solutions satisfy 

IK.fc - u\\L°°({(x,t): \d(x,t)\>VZ}) = C(e k+2 ), as e -> 0, 
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with u the sharp interface limit defined in (12. U via the volume preserving mean curvature 
flow (11.12b starting from Tq. We can therefore interpret Theorem 12. 3 l as a result of conver- 
gence of the mass conserving Allen-Cahn equation ( II . lb to the volume preserving mean 
curvature flow (11.12) : 

sup ||u e (-,t)-fi(-,t)||^ ( n) = 0(e 1/4 ), ase^O. 

0<t<T 

Organization of the paper. The organization of this paper is as follows. In Section [3] we 
present the needed tools which are by now rather classical. In Section|4] we perform for- 
mal asymptotic expansions of the solutions (u e (x, t), X e (t)). This will enable to construct 
the approximate solutions (u e ,k(x, t), \ E ,k(t)), and so to prove Theorem l2.21 in Section|5] 
Last we prove the error estimate of Theorem 12.3 l in Section |6] In particular and as men- 
tioned before, a precise understanding of the error between the actual and the approximate 
Lagrange multipliers will be necessary (see subsection l6.U . 

Remark 2.4. Through the paper, the notation ip £ sa J2i>o e ^ represents asymptotic ex- 
pansions as e — > and means that, for all integer k, ij) e = Xa=o + 0{£ k+1 ). 

3. Preliminaries 

For the present work to be self-contained, we recall here a few properties which are 
classical in the works mentioned in the introduction, ]25], 0], fl24), 02), (Til. 031. flUl. 
and the references therein. 

3.1. Some related linearized operators. We denote by 6o(p) := tanh(^=) the standing 
wave solution of 

9 Q " + f(e Q ) = onl, 
o (-oo) = -l, O (O) = O, fl (oo) = l, 

which we expect to describe the transition layers of solutions u e observed in the stretched 
variable. Note that, for all m £ N, 

(3.1) D™[6 (p) - (±1)) =<D(e-^) asp^ioo. 

We then consider the one-dimensional underlying linearized operator around 9q, acting 
on functions depending on the variable p by 

(3.2) Cu := -u pp - f'(6 (p))u. 

Lemma 3.1 (Solvability condition and decay at infinity). Let A(p, s, t) be a smooth and 
bounded function on K x U x [0,T], with U C R w_1 a compact set. Then, for given 
(s, t) e U x [0, T], the problem 

C^j := -tP pp - f(6 (p))ilJ = A(p, s,t) on E, 

has a solution (which is then unique) if and only if 

(3.3) / A(p,s,t)6 '(p)dp = 0. 

JR 

Under the condition ( 13.31 ), assume moreover that there are real constants and an 
integer i such that, for all integers m, n, I, 

(3.4) D™D?D l t [A(p,a,t)-A ± ] = 0(\p\ i e- y/5M ) asp^±oo, 
uniformly in (s,t) G U x [0, T]. Then 

(3.5) D?D?D[[1,(p,8,t)-j^]=OQ P \ i e-^>'\) asp -+ ±oo, 
uniformly in (s,t) G U x [0,T]. 
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Proof. The lemma is rather standard (see J4], among others) and we only give an outline 
of the proof. Multiplying the equation by 6q and integrating it by parts, we easily see that 
the condition d3.31 l is necessary. Conversely, suppose that this condition is satisfied. Then, 
since 6q is a bounded positive solution to the homogeneous equation ip pp + f'(9o(p))tp = 
0, one can use the method of variation of constants to find the above solution tp explicitly: 

TP(p,s,t) = -6 '(p) f (e Q '- 2 (0 [°° A(Z,8,t)B '(p)dt) dC 



JO \ J( / 

Using this expression along with the estimates (13.4b and (13.11 ). one then proves ( 13.5b . □ 

Note also, that after the construction of the approximate solutions u E: k, we shall need the 
estimate of the lower bound of the spectrum of a perturbation of the self-adjoint operator 
—A — e~ 2 /' (u e .k) proved in [131 . This will be stated in Section[6] 

3.2. Geometrical preliminaries. The following geometrical preliminaries are borrowed 
from 0141 . to which we refer for more details and proofs. 

Parametrization around T. As mentioned before, we call T = Uo< t <T(^t x {£}) me 
smooth solution of the volume preserving mean curvature flow ( 11.12b . starting from To; 
we also denote by il t the region enclosed by T t . Let d be the signed distance function to T 
defined by 



(3.6) d(x,t) 



-dist(a;,r t ) for a; € fit 
dist(x, T t ) for x <G ft \ Tk. 



We remark that d is smooth in a tubular neighborhood of V, say in 

M 3S (T t ) := {x e fi : \d(x,t)\<35}, 

for some S > 0. We choose a parametrization of T t by Xq(s, t), with s S U C R^ -1 . We 
denote by n(s, t) the unit outer normal vector on dVl t = IV For any < t < T, one can 
then define a diffeomorphism from (—3(5, 3d) x U onto the tubular neighborhood A3i(r t ) 
by 

X(r,s,t) = X (s,t) + rn{s,t) = i£ Af 3 s(T t ), 

whose inverse is denoted by r = d(x,t), s = S(x,t) := (S 1 (x, t), ■ ■ ■ , S N ^ 1 (x,t)). 
Then Vd is constant along the normal lines to T t , and the projection S(x,t) from x on 
T t is given by Xq(S(x, t),t) = x — d(x, t)Vd(x, t). For x — Xo(s, t) £ T t denote by 
Ki(s,t) the principal curvatures of r t at point a; and by V(s,t) := (Xo) t (s,t).n(s,t) the 
normal velocity of T t at point x. Then, one can see that 

JV-l 

(3.7) K(s,t) := K »( s ^) = &d(X (s,t),t), 

i=l 



N-l 

(3.8) h{s,t) :=-J2 K * 2 M) = -{Vd.VAd){X {s,t),t), 

i=l 

(3.9) V(s,t) := (X ) t (s,t).n(s,t) = -dt{X{r,s,t),t). 

In particular, dt{x,t) is independent of r = d(x,t) in a small enough tubular neighbor- 
hood of r t . Changing coordinates form (x, t) to (r, s, t), to any function <f>(x, t) one can 
associate the function <j>(r, s, t) by 

4>(r,s,t) = 4>(X (s,t) + rn(s,t),t) or <j>(x,t) = <f>{d{x, t), S(x, t),t). 
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The stretched variable. In order to describe the sharp transition layers of the solutions u £ 
around the limit interface, we now introduce a stretched variable. Let us consider a graph 
over r t of the form 

Tf = {X(r,s,t) : r = eh £ (s,t),se U}, 

which is expected to represent the level set, at time t, of the solutions u e . We define the 
stretched variable p(x, t) as "the distance from x to Tf in the normal direction, divided by 
e", namely 

/ A d(x,t)-eh £ (S(x,t),t) 

(3.10) p[x,t):= . 

£ 

In the sequel, we use (p, s, t) as independent variables for the inner expansions. The link 
between the old and the new variable is 

x = X(p, s, t) := X (e(p + h £ (s, t)), s, t) = X (s, t) + e(p + h £ (s, t))n(s, t). 

Changing coordinates form (x, t) to (p, s, t), to any function ip(x, t) one can associate the 
function i/j(p,s, t) by 

(3.11) $(p,*,t) =ii(Xo(a,t)+e{p + h e {a,t))n{a,t),t), 

or ij)(x, t) = $( *lB&=2MSi*&£ : s(x, t),t). A computation then yields 
£ 2 {d t ^ -Atl>) = - 4> pp - e(V + Ad)^ 

(3.12) + e 2 [df^ - A r 7/> - {dfh £ - A F h £ )ijj p } 

+ e 2 [2V r /i £ .V r ^- \W r h s \% p }. 

where 

JV-l JV-l JV-l N-1 

Of : », ■ Y. Std s i, V r := ^ S ' n * ■ AF : E AS ' iK + E V^.V^. 

i—l i—1 i—1 i-.j — l 

Here Ad is evaluated at (x, t) = (X (s, t) + e(p + h e (s, t))n(s, t), t), so that (13.7b and 
( 13.8b imply 

Ad = Ad{X (s, t) + e{p + h s (s, t))n(s, t),t) 
(3 - 13) ~ n(a, t) - s{p + h £ {s, t))h(a, ~ E + t))%(a, t), 

i>2 

where 6j(s, t) (i > 2) are some given functions only depending on T t . 
Last, define 

eJ%p,s,t) :=dX(p, S ,t)/d(p 7 s) 

the Jacobian of the transformation X so that, in particular, dx = eJ e (p, s, t) dsdp. Then, 
for all p e R, s e U and < t < T, we have 

JV-l 

(3.14) J S (p,S,t)= l[[l + E(p + h B (8,t))Ki(8,t)]. 

i=l 

4. Formal asymptotic expansions 

In this section, we perform formal expansions for the solutions u £ (x,t) of ( 12.3b . We 
start by outer expansions to represent the solutions "far from the limit interface", then make 
inner expansions to describe the sharp transition layers. Last, expansions of the nonlocal 
term X £ (t) are performed. In the meanwhile we shall also discover the expansions of the 
correction terms h £ (s,t) defined in ( 13.10b . 

We assume that the solutions u £ (x, t) are of the form 

(4.1) u £ (x, t) sa uf(t) := ±1 + euf(t) + e 2 u J (<) + •• • (outer expansions), 
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for x G fit (corresponding to x G fi \ fit (corresponding to uf(t)), and away 

from the interface r t , say in the region where \d(x, t)\ > \fe as we expect the width of 
the transition layers to be 0(e). Near the interface T t , i.e. in the region where \d(x, t)\ < 
y/e, we assume that the function u e (p, s, t) — associated with u e (x, t) via the change of 
variables (13. 1U — is written as 

(4.2) u 6 (p, s, t) ss Uo(p, s, t) + eu\{p, s, t) + e 2 U2(p, s,t) + ■■ ■ (inner expansions). 

We also require the matching conditions between outer and inner expansions, that is, for 
all i G N, 

(4.3) Ui(±oo, s,t) = uf(t) (matching conditions), 

for all (s, t) G U x [0, T]. As we expect the set p = to be the level set of the solutions 
(see subsection !3.21 i we impose, for all i G N, 

(4.4) Ui(0, s,t)=0 (normalization conditions), 

for all (s,t) G U x [0,T]. 

As far as the nonlocal term A e (t) is concerned we assume the expansions 

(4.5) X s {t) w X (t) + eAi(t) + e 2 \ 2 {t) H (nonlocal term). 

Last, the distance correcting term h e (s, t) is assumed to be described by 

(4.6) eh E (s, t) m eh%(s, t) + e 2 /i 2 (s, £) + ■■• (distance correction term), 

for all (s,t) G U x [0,T]. 

In the following, by the (complete) expansion at order 1 we mean 

{d(x,t), Ao(t),Ui(/J, s,t),u 1 (t)} (expansion at order 1), 

and by the (complete) expansion at order i > 2 we mean 

(4.7) {hi-i(s, t), Xi-x(t), Ui(p, s, t),uf(t)} (expansion at order i > 2). 
Let us also recall that we have chosen 

f(u)=u(l-u 2 ), W(u) = ±(l-u 2 ) 2 . 

4. 1 . Outer expansions. By plugging the outer expansions (14. U and the expansion (14.51 l 
into the nonlocal partial differential equation (12. 3K we get 

(4.8) = - (4 (t)) 3 - e\ e (t)(l ~ (uf(t)) 2 ). 

Since uf(t) w 2~)i>o elu t(t)> wnere w o(^) = ±1> an elementary computation yields 

-eA e (t)(l- (u±(t)) 2 ) f E A pW E 

and 

kW«e( E «p(*) E u^wufc*))^- 

i>0 yp+<j=i k+l=q J 

Hence, collecting the e terms in ( 14.8l l. we discover = (t) — 3uf(t)(v,Q (t)) 2 so that 
= 0. Next, an induction easily shows that 

uf(t) = foraUi>l. 

Therefore the outer expansions are already completely known and are trivial: 

(4.9) uf(t)=±l. 

In other words, thanks to the adequate form of the Lagrange multiplier, the outer expan- 
sions are independent of the expansion of the nonlocal term. This is in contrast with the 
equation considered in |[T4l . 
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4.2. Inner expansions. It follows from (13.12b that, in the new variables, equation (12.31 ) is 
recast as 

(4.10) u epp + u £ - {u £ f = e\ e (t){l - {u e f) - e(V + Ad)u ep 

+e 2 [dfu E - A r u E - (dfh E - A T h E )u Ep ] 
+e 2 [2W r h E .W r u Ep - \V r h £ \ 2 u epp ]. 

The e° terms. By collecting the e° terms above and using the normalization and matching 
conditions (14.3b . ( 14.4b we discover that uq(p, s,t) = 0q{p), with 8q the standing wave 
solution of 



(4.11) 



0o" + /(0o) = O onR, 

o (-oc) = -l, O (O) = O, o (oo) = l. 



Formally, this solution represents the first approximation of the profile of the transition 
layers around the interface observed in the stretched coordinates. Note that since /(it) = 



u — u 3 , one can even compute #o(p) = tanh(^). 



The e 1 terms. Next, since u £ (p, s, t) as J2i>o u i(P' s > O 6 *' wnere uo(p, s, t) = 0o(p), an 
elementary computation yields 

(4.12) z\ s (t)(l-(u e f(p,s,t)) I X P (t)(3 q (p,s,t)\e t+1 , 

i>0 \p+q—i J 



where 

2 



Pq(j>,S,t) 

and also 



o if? = 



(4.13) {u E ) 3 {p,s,t) as ^ u p (p,s,t) ^ u k (p,s,t)ui(p,s,t)\e\ 

i>0 \p+q=i k+l=q J 

Hence, plugging the expansion ( 13.131 ) of Ad into ( 14.10b and collecting the e terms, we 
discover 

(4.14) Cm := -u lpp - f'(6 (p)) Ul = (V + «)(«, t)0 o '(/>) - (1 - # 2 (p))A (i). 

For the above equation to be solvable (see Lemma [37X1 for details) it is necessary that, for 
all (s,t) e U x [0,T], 

/ Cui(p,s,t)6 '{p)dp = 0, 

which in turn yields 

?0 )O0 



(4.15) V(s, t) = -k(s, t) + aX (t), 



As seen in subsection !3.2l the above equation can be recast as 

(4.16) d t (x,t) = Ad(x,t) - crAo(t) for x 6 T t . 

Now, in view of gJU), we can write = X! 00 (^o" + /(^o))6 , o' = fl^o" ~ W'{6 Q ))0 ' 
and find the relation 1 - 6 2 = V2do, so that a = \f2. Plugging this and ( 14.15b into ( 14.14b 
we see that Cu\ = 0. Therefore, the normalization Mi(0, s, t) = implies 

(4.17) Ul {p,s,t) = 0. 
Again this is in contrast with the equation considered in 0141 . 
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The e % terms (i > 2). Now, taking advantage of Uo(p, s,t) = (p) and of Ui(p, s, t) = 
we identify, for i > 2, the e l terms in all terms appearing in (14.101 . In the sequel we omit 
the arguments of most of the functions and, by convention, the sum is nu U if b < a. 
Using ( 14.131 ) we see that the e % term in u epp + u e — (u £ ) 3 is 

i-2 i-2 

(4.18) - Cui - 9 ^2 u k Ui- k - ^ u p UkUl < - term 

k=2 p=2 k+l=i-p 

In view of ( 14.121 ), the e l term in e\ £ (t)(l — (u e ) 2 ) is 

(4.19) Ai-ifl-flo 2 )- J2 X p J2 UkUl (term2). 

p+q— 2— l,g^0 k+l—q 

In order to deal with the term — e(V + Ad)u £p , we first note that ( 13.13b and ( 14.6b yield 
the following expansion of the Laplacian 

(4.20) Ad^ n-J2(b 1 h i + Si) e\ 

i>i 

with 

i 

(4.21) 5 i = 5 i (p,s,t) = Y,Ck(s,t)p k 

k=0 

a polynomial function in p of degree lower than i, whose coefficients c/;(s,t) are them- 
selves polynomial in (hi , hi-i) which are part of the formal expansions at lower orders, 
and in (&i, ...&,) which are given functions. Among others, we have 5\(p, s,t) = &i(s, t)p 
and 5 2 {p, s, t) = 6 2 (s, t)(p + hi(s, t)) 2 . Combining u ep « 8 ' + e 2 u 2p + • • • and ( 14.201 ), 
we next discover that the e l term in — e(V + Ad)u ep is 

(4.22) hhi-tOo' + Si-rfo' - (V + k) u (i _ 1)p 

+ Ep=i ( b i h P + S p) u (i-i~ P )p ( term 3). 

We see that the e l term in e 2 [dfu £ — A r u e — (dfh e — A T h e )u ep ] is given by 

i-3 

(4.23) (df - A T ) Ui - 2 - (df - A I > J _ 1 # ' - J2(df - A T )h p u (l ^ p)p (term 4). 

P =i 

Note that 

\eV T h E \ 2 « e 2 |V r /n| 2 + (2V r /n.V r ^_! + m ) e\ 

i>3 

where 

rji = r)i(s,t) := ^ V r /i p+ i(s,i).V r /i g+ i(s,t) 

p+<3=i-2 : p^O : q#0 

depends only on the derivatives of hi,...,hi-2- Combining this with u epp rj Oq" + e 2 u 2pp + 
■ ■ • , we discover that the e l term in — e 2 \ V T h e \ 2 u Epp is 

i-3 

(4.24) - &(V r /ii.V%_i)(9 " - |V r /n| 2 u (j _ 2)pp - a kUkpp (term 5), 

k=0 

where otk = otk{s, t) depends only on the derivatives of h\,...,hi-2 and j3 2 = 0, fii = 2 if 
i > 3. 

Last, since V r w ep ss e 2 W T U2 P + ■ • • , we see that the e l term in e 2 [2V r h £ .'V r u ep ] is 

i-2 

(4.25) 2^ V r /i l _i_ fc .V r u fcp (term 6). 

k=2 
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Hence, in view of the six terms appearing in (14.181 . (14.19) , (14.22b . ( 14.23b . ( 14.24b . ( 14.25b . 
when we collect the e l term (i > 2) in d4.10b we face up to 
(4.26) 

Cm = (M r h^ 1 )6 Q l -(l-9 a 2 )\^ 1 +l3 l (V r h 1 .V r h^ 1 )9 a // + \V r h 1 \\ { ^ 2)pp + R l ^ 
where A4 r denotes the linear operator acting on functions h(s, t) by 

(4.27) M r h := dfh - A r h - hh, 

and where Ri-i = Ri-i(p, s,t) contains all the remaining terms. Since it is important 
that Ri-i does not "contain" we have to leave |V r /ii \ 2 U( i _ 2 ) PP for the case i = 2, 
but with a slight abuse of notation we can "insert" |V r /ii| 2 "(;-2)pp in R4-1 for i > 3. As 
an example, for i = 2 we see that 

(4.28) Rx{p,s,t) = -8^,3,1)60' \p) = -&i(s,t)p<V(p), 
so that we infer that, for all integers m, n, I, 

(4.29) D™ 1 D l t [Ri{p, s,t)] = 0(\p\e~^ p \) as p -> ±00, 

uniformly in (s, t). Now, for i > 3, we isolate the "worst terms" — which are the Si's — in 
Ri-i and write 

i-3 

(4.30) = -5i-i6 ' - ^2 M(»-i-p)p + r,_i, 
where r,_x = r.;_i (p, s, t) contains all the remaining terms. 

Lemma 4.1 (Decay of i?i-i). Lef i > 2. Assume that, for any 1 < k < i — 1, there holds 
that, for all integers m, n, I, 

(4.31) C™^^K(p, S ,t)] = 0(|p| fc - 1 e -^l' , l) asp^±^, 
uniformly in (s,t) £ U X [0,T]. Then, for all integers m, n, I 

(4.32) D-D«D|[i? i _ 1 (p,s,t)]=0(|pr- 1 e-^l' , l) as p -> ±00, 
uniformly in (s,t) € E/ x [0, T]. 

Proof. Let us have a look at expression (14.30b of J2,_i. By a tedious but straightforward 
examination we see that r^_i (p, s, t) depends only on 

• V(s, t), k(s, t), bi(s, t), bi(s, t) which are bounded given functions 

• A (t),...,A i _a(t) 

• hi(s, t), hi_2(s, t) and their derivatives w.r.t. s and t 

• uo(p 1 s 1 t) = 8(p),ui(p, s,t) = 0, Ui-i(p, s, t) and their derivatives w.r.t. p, 
s and t 

in such a way that it is 0{\p\ l ~ 2 e~^ 2<[P<[ ) as p — s- ±00. Concerning the term 

i-3 

-5i-x{p, S, t)9 '(p) - 2J MA S > W (i-l-p)p: 

p=l 

the fact that it behaves like d4.32b follows from ( 14.31b and the fact that S p (p, s, t) grows 
like \ p\ p , as seen in d4~2Ti □ 

Now, in virtue of Lemma [3~Tl the solvability condition for equation ( 14.26b yields, for all 

(M), 

(M r h i . 1 )(s,t) f 9 ' 2 -h-i(t) f (1-9 Q 2 )9 Q '+ ( Ri-i(;8,t)6 ' = 0. 



Note that the term — /3j(V r /ii.V r /ij_i)#o" does not appear above since L6q"6o' = 0. 
Note also that the term — \V T hi\ 2 U( i _ 2 ) pp does not appear for the same reason if i = 2, 
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and because it can be "hidden" in R4-1 for i > 3 without altering the fact that R4-1 does 
not depend on hi—i. The above equality can be recast as 

(4.33) {M T h i - 1 )(s,t) = aX i - 1 (t)-a* [ R l - 1 (p,s,t)0 o / {p)dp, 

with a defined in ( 14.15b and a* := (J R 6 ' 2 )~ . Note that, thanks to 1 - 9 2 = y/29 ', we 
have a = \/2 (as seen before) and also a* = |vZ 

Let us have a look at i = 2. From (14.28b and the fact that j ' p^' (p) dp — (odd 
function), we see that ( 14.33b reduces to 

(4.34) (M r h 1 )(s,t) =aXx(t). 

Assume that hi satisfies the above equation. Then since u\ = trivially satisfies ( 14.31b . 
Lemma l4~T1 implies that R\(p, s, t) together with its derivatives are C ( | /o| e ^/^'l' 3 ! ) as p — >• 
±00. It follows from Lemma lXTI fhat 

(4.35) Cu 2 = (M r h 1 )9 Q l - (1 - 9 2 )X 1 + |V%| 2 O " + Ri, 

admits a unique solution U2(p, s, t) such that 1*2(0, s, t) — 0, which additionally satisfies 
D™D?D l t [u 2 (p,s,t)] = 0(\p\e-^). 

Now, an induction argument straightforwardly concludes the construction of the inner 
expansions. 

Lemma 4.2 (Construction by induction). Let i > 2. Assume that, for all 1 < k < i — 1 

the term Uk is constructed such that 

(4.36) D™D^D l t [u k (p,s 1 t)}=0(\p\ k - 1 e-^) asp^±^, 

uniformly in (s,t) G U x [0,T]. Assume moreover that hi-i(s,t) satisfies the solv- 
ability condition ( 14.33b . Then one can construct Ui(p,s,t) solution of ( 14.26b such that 
Ui(0, s, t) = and 

(4.37) D™D n 8 D l t [u i { P ,s,t)]=0{\p\ i - 1 e-^) as p -> ±00, 
uniformly in (s, t) € U X [0, T]. 

Remark 4.3. Note that the cancellation Hi = implies that the term Ui(p, s, t) appearing 
in the expansion of the solutions of ( II. lb behaves like (D(\p\ l ~ 1 e~^ 2 ~\ p \), where the term 
Ui(p, s, t) appearing in the expansion of the solutions of (II. lib behaves like 0{\p\ l e^^ p ^) 
(see Q3|). 

4.3. Expansions of the nonlocal term X e (t) and the distance correction term h e (s, t). 

By following lfT4l subsection 5.4] with y^i playing the role of 8, we see that an asymptotic 
expansion of the conservation of the mass (11.5b yields 

(4.38) = — / u e (x, t) dt^h + I 2 + I 3 , 

dt J n 

where I\ = 0, since in our case uf (t) = ±1, and 

(4.39) h ■= I dlu E {p,s,t)e.J £ (p,s,t)dpds, 

J\p\<l/Ve 

(4.40) h ~ f {-V -edfh E ){s,t)d p u E {p,s,t)J e {p,s,t)dpds 7 

J\p\<l/Ve 

Combining df := d t + J^iLi 1 Std s i with Uo(p, s, t) = (p) and u\(p, s, t) = 0, we 
see that 

N-l 

dfu E (p,s,t) w ^e l [d t + Sfd a k]ui(p,s,t). 

i>2 k=l 
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In view of the above inner expansions, this implies 

^tf e (p,8,t)«53 e<0 (H < " le "^ |p| )' 



i>2 



where O 



into ( 14.391 1, we get 



(ipl^ 



-V2\p\ 



depends only on expansions at orders < i — By plugging this 



i>3 



where 7;_2 = (t) depends only on expansions at orders < i — 2. 
We now turn to the term 1%. We expand 

(-V-e$h e ) (s,t) ps d t (X (s,t),t) -^V^/nfot), 



i>i 



and 



d p u e (p,s,t) ps o '(p) +yVc? p i<t(p i s,*)- 

Expanding the Jacobian ( 13.14b and using ( 13.71 ), we get 

J e (p, s, t) ps 1 + Ad(X (s, t), f) e(p + /i £ (s, t)) + e'lH-i, 



i>2 



where fc-i = fM-i (p, s, t) depends only on expansions at orders < i — 1. Multiplying the 
three above equalities, we see that the integrand in J3 expands as 

60 d t + e6 ' [-df^ + Ad + pd t Ad] + ^ £ 4 (9 '(-<9 t r /i 4 + MtAd + 

i>2 

where = Uj_i(p, s,t) depends only on expansions at orders < i — 1. We inte- 

grate this over s £ U and |p| < l/y/e and, using S\ p \<.x/^/E^ ~ Jr^o'(p)^P = 2 and 
/|p|<i/V£ P@o (p) dp — (odd function), we discover 



d t (s,t)ds + e / (-dfh 1 + (d t Ad)h 1 )(s,t)ds 
u 



(-dfhi + (d t Ad)hi)(s,t) ds + w;_i 



where = depends only on expansions at orders < i — 1, Using (14.161 1 to 

substitute d t , ( 14.34l i to substitute dfhi, ( 14.33b to substitute we have 

i/ 3 ps / (Ad - ctA ) ds + e [ (-A r /ii - bi/j x - ctAi + (d 4 Ad)/ii) ds 
2 7 c; J(7 



E< 

j>2 



(-A r /i i - 61^ - uAj + (d t Ad)hi) ds + Ci-i 



where = £j_i(t) depends only on expansions at orders < i — 1. 
Last, using Jy A r /l, dS = 0, we see that I2 + I3 ~ reduces to 



(4.41) 
(4.42) 
(4.43) 



aXo{t) = Ad(;t) 



aAi(t) = -[&i(-,i)-d t (-,t)Ad(-,i)]M-,*) 



<rAi(t) = -p»i(-,t) -d t (;t)Ad(;t)]hi(;t) + A*_i(t) (» > 2), 



where </>(•) := Tyr J v 4> denotes the average of <j> over T t (parametrized by U), and A{_i (t) 
depends only on expansions at orders < i — 1. Moreover if we plug (14. 4U , (14.421 and 
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( 14.43l > into ( 14.16b . (14.34b and ( 14.33b . we have the following closed system for d, hi,.., hi 

onU x [0,T]: 



(4.44) dt = Ad - Ad[; t) 

(4.45) dfh x = A r h ± + bih t - [&!(•,*) -dt(-,*)Ad(-,t)]/ii(-,t) 



(4.46) %ht = A hi + bi hi - [h(;t) - d t (;t)Ad(;t)]hi(;t) + A^_ x (t) (i > 2). 

5. The approximate solutions u E) fc, A £ .fc 

In order to construct our desired approximate solutions and prove Theorem 12.21 let 
us first explain how the previous section enables to determine, at any order, the outer 
expansion ( 14.1b . the inner expansion (14.2b . the expansion of the nonlocal term (14.5b . and 
the expansion of the distance correction term (14.6b . 

First, as seen before, the outer expansion (14.1b is already completely known since 
uf(t) = OforalH > 1. 

Recall that T = Uo<t<T(rt x {t}) denotes the unique smooth evolution of the volume 
preserving mean curvature flow (11.12b starting from To CC O, to which we associate the 
signed distance function d(x, t). Hence, defining Xo(t) as in (14.41b and u\(p, s, t) = as 
in ( 14.17b . we are equipped with the first order expansion 

(5.1) {d(x,t),\ (t),ui(p,s,t) = 0}. 

Next, since T t is a smooth hypersurface without boundary, there is a unique smooth so- 
lution hi(s,t) to the parabolic equation ( 14.45b . Assuming hi(s, 0) = for s 6 U, we 
see that hi(s,t) = 0, which combined with (14.42b yields Xi(t) = 0. Notice that these 
cancellations are consistent with the observation of [8 1 that " (II. lb has better volume pre- 
serving properties than the traditional mass conserving Allen-Cahn equation ( 11.1 lb ". In 
Section |4] we have defined U2(p,s,t) as the solution of (14.35b . which now reduces to 
£i*2 = — &i(s, t)p9o'(p). This completes the second order expansion, namely 

(5.2) {MM) = 0,Ai(t) = 0,u 2 {p,s,t)}. 

Now, for i > 2, let us assume that expansions {hk-i(s, i), Xk-i(t),Uk(p, s, t)} are con- 
structed for all 2 < k < i. Therefore we can construct Aj_x(t) appearing in ( 14.46b . 
Assuming hi(s, 0) = for s E U, there is a unique smooth solution hi(s, t) to the para- 
bolic equation (14.46b . This enables to construct Xi(t) via ( 14.43b . Now, hi(s, t) satisfies the 
solvability condition ( 14.33b at rank i, so that Lemma l4~2l provides Ui+x(p, s, t), the solu- 
tion of (14.26b at rank i + 1 with m+i(0, s, t) = 0. This completes the construction of the 
i + 1-th order expansion {/ij(s, t), Xi(t),Ui + i(p, s,t)}. 

Note also that, from the above induction argument, we also deduce the behavior (14.37b 
for all the s, t)'s. 

Proof of Theorem 12.21 We are now in the position to construct the approximate solutions 
as stated in Theorem l2.2l Let us fix an integer k > max(N, 4). We define 



ps,k(x,t) 



k+2 

d(x,t) - e l hi{S(x,t),t) 



i=l 



£ 



k+3 



u%(x,t) := e {p E . k (x,t)) +^e i u l {p^ k {x,t),S{x,t),t), 

»=i 

uZ%{x,t) :=u(x,t), 

k+2 



X £ , k (t) :=X (t) + 
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where u is the sharp interface limit defined in (12.1) . We introduce a smooth cut-off function 
C(z) = (s(z) such that 

'C(z) = l X\z\<y/i, 

C(z) = if \z\ > 2y/E, 

0<z('(z)<4 if y/e < \z\ < 

For x <S f2 and < t < T, we define 

:= C(dOM)K? fe CM) + [1 - ((d(x,t))]u™£(x,t). 

If e > is small enough then the signed distance d(x, t) is smooth in the tubular neighbor- 
hood -A/" 3% /j(r), and so is u™ k (x, t). This shows that u* k is smooth. 

Plugging (it* k (x, t), A e (t)) into the left hand side member of (12.5b . we find a error term 
6* k (x, t) which is such that 

• S* k {x 7 t) = 0on{\d(x,t)\ > 2y/e} since, then, u* fc = = ±1, 

• ll^fclU 00 = C(e fc+2 ) on {\d{x,t)\ < ^/e} since, then, u* e k = u™ k and the 
expansions of Section|4]were done on this purpose, 

• ll^fclU 00 — 0{e k ' ), for any integer k* , on {y/e < \d(x,t)\ < 2-y/i} since, then, 
the decaying estimates ( 13.11 ) and ( 14.37b imply that u* k — u° u k = u* k — ±1 = 

0(e 2 v^ ), valid also after any differentiation. 

Hence \\8* fc||_L°°(six(o.T)) = 0(e k+2 ), which is even better than ( 12.5b . Also u* k clearly 
satisfies ( 12.7b . 

Now, to ensure the conservation of the mass of the approximate solutions, we add a 
correcting term (which depends only on time) and define 

u £ ,k{x,t) := u* e k (x,t) + p| J^(u* k (x,0) - u* tk (x,t)) dx, 

which then satisfies (12.8b . and still d2.7b . Note also that subsection 14.31 implies that the 
correcting term 



0) - u* Jx,t))dx = - 



/ / dtU* k (x,T) drdx 

JQ JO 



is 0(e k+2 ) together with its time derivative. Hence, when we plug u e> k = u* . + 0(e k+2 ) 
into the left hand side member of (12.5b . we find a error term S £t k whose L°° norm is 

o(e k ). a 



6. Error estimate 

We shall here prove the error estimate, namely Theorem 12.31 For ease of notation, 
we drop most of the subscripts e and write u, A, Uk, A^, 8k for u e , A e , u e ^, A £i /-, 5 £) k 
respectively. By || • ||, || • |2+ p we always mean || • ||x, 2 (0) 5 II ' \\l 2 +p(Q) respectively. In 
the sequel, we denote by C various positive constants which may change from places to 
places and are independent on e > 0. 

Let us define the error 

R(x, t) := u(x, t) — Uk(x, t). 

Clearly ||i£||i,oo < 3. It follows from the mass conservation properties (11. 5) , (12.8b . and the 
initial conditions ( 12.10b that 

(6.1) [ R(x,t)dx = Q forallO<<<T, \\R(-, 0)|| = 0(e k ^^). 
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We successively subtract the approximate equation (12.5b from equation Jl.lt . multiply 
by R and then integrate over ft. This yields 



(6.2) 

where 
(6.3) 



+ \ [ (f(u)-f(u k )-f(u k )R)R- f S k R-±A, 
£ Jn Jn £ 

A = A(t) := / [eA(l -u 2 )R-e\ k {l - u k 2 )R]. 
Jn 



Since (f(u)-f{u k )-f'{u k )R)R = -3u k R 3 -R 4 = 0{R 2+P ), where p := min(A 1), 
we have 



4 / (f(u)-f(u k )-f(u k )R)R 



^^cwRwixi^^c.wRinvRw 2 , 

where we have used the interpolation result lfT4l Lemma 1]. We also have \J n 5 k R\ < 
||4||co||i?|| =0(e fe )||i?||,sothat 



(6.4) 



\R\A\R\\<- I \VR\ 2 + ± I f{u k )R 2 
«£ Jn £ Jo 



l-dWRinVRf + 0(e k )\\R\\ - -1a. 



We shall estimate A in the following subsection. As mentioned before, this term is the 
main difference with the case of a strictly nonlocal Lagrange multiplier: its analogous for 
equation ( 11.1 U is (eX — eX k ) J n R which vanishes, see Ifl4l . 

Since k > max(N, 4) we have k — \ > - = . , 4 4 1N , so that the second estimate in 

v ' ' 2 p mm(-^,l)' 

( 16. U allows to define i e > by 

(6.5) t e :=sup{i >0,V0<r <t, \\R(-,t)\\ < {2C 1 y l/p e i/p } . 

We need to prove that t e = T and that the estimate 0(e 4 ^ p ) is actually improved to 
Oie k ~ 2 ). In the sequel we work on the time interval [0, t e ]. 

6. 1 . Error estimates between the nonlocal/local Lagrange multipliers. It follows from 
( 12. 9i that the term A under consideration is recast as 

(6.6) A= A E _A± Ek+ 0(^) E ^ 

ts ts k ts k 

where 

A k = A k (t) := f /(u fc ), B k = B k (t) := [ 1 - u k 2 , E k = E k (t) := [ (1 - u k 2 )R, 
Jn Jn Jn 

and A, B, E the same quantities with u in place of u k . 
Lemma 6.1 (Some expansions). We have, as e — >• 0, 

A k =e 2 a + 0(e 3 ), B k = ep + 0(e 2 ), 

where 

N-l 

a = a(t) := / V Ki(s,t)ds / pf(8 (p))dp, (3 := 2y/2\U\, 
Ju i=1 Jr 

and 

E k = OWi\\R\\). 
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Proof. We have seen in Section[5]that Uk = lit +0(e k+2 ) so it is enough to deal with A* k 
B k and El*. The lemma is then rather clear from the expansions of Section|4] We have 

A* k = [ f(u* k )(x,t)dx= [ f(u* k )(x,t)dx + 0(e-^) 



|d(x,t)|<2v/e J\d(x,t)\<^e 

= 11 f(9o(p) + 0(e 2 ))eJ £ (p, S ,t)dsdp + 0(e-^). 

Ju J\p\<l/Vs 

Using .P(p lSl t) = 1 + epEtT 1 + ^(e 2 ) and f\ P \<i/^ fiW) dp = (odd 

function), one obtains the estimate for AT. The estimate for B k follows the same lines and 
is omitted. Last, the Holder inequality yields \E k \ < (J n (l-u k 2 ) 2 ) 1/2 \\R\\ = 0{^e\\R\\) 
since, again, dx = eJ £ (p, s, t) ds dp. □ 



As a first consequence of the above lemma, it follows from (16. 6t that 

(6.7) k= ± E -^E k +0{e k+ H\R\\. 

t> a k 

Next, in view of the above lemma, u = u k + R and \\R\\ = C(e 4 / p ), we can thus perform 
the following expansions 



since Jn R = o. 



B- 1 = B k 



A = A k + [ (1 - 3u k 2 )R - 3 f u k R 2 - I 
Jn Jn Jn 

= A k + 3E k -3 [ u k R 2 + 0{\\R\\lX P p 
Jn 

B k B k 



R 3 

n 



kti. \ 2 

B k B k V B k 



Bk -i L + 2j^g , In* 2 , p/a^ ^ + Q ( \\ R \\ 



and 



E = E k -2 [ u k R 2 + 0(\\R\\1%1). 
Jn 

It follows that, using E k = 0(*/e\\R\\) and A k = 0(e 2 ) (see LemmaloTTTi, 

AE=A k E k -2A k f u k R 2 + 0(e 2 \\R\\ 2 2 l P p ) + 3E k 2 + 0(V^\\Rf) 
Jn 

+ 0(^\\R\\ \\R\\t$ + 0(V^\\R\\ 3 ) + 0(||i?|| 4 ) + C{\\R\\ 2 \\R\\l + + D 
+ 0(\\R\\lX P p V~e\\R\\) + 0(\\R\\lX P p W R W 2 ) + °(W R WtXp) 

=A k E k + 3E k 2 -3E k f u k R 2 -2A k f u k R 2 + 0(e 2 \\Rf 2 X P p ) + 0(\\R\\ 3 ), 
Jn Jn 

since ||i?||2+p = C(ll-^l| 2 )- Then, using Lemma loTTl and a methodical grading of the C's 
(the typical used arguments being \\R\\lXp = 0{\\R\\ 2 ) and, e.g., ||i?j| 4 = 0(^/e\\R\\ 3 ) by 
the definition of t e ), we arrive at 

± E - ^E k =B k - 1 
B B k 



3E k z - 3E k I u k R z + -^-E k I 2u k R 

2+P\ 



2A k / u k R 2 + 0(e 2 \\R\\tX P P )+0(\\R\\ 3 ) 



IS 



M. ALFARO AND P. ALIFRANGIS 



which in turn implies 



A A k 
—E -E k 

B B k 



3E k 2 - 3E k J Q u k R 2 + j±E k J n 2u k R 
B k 

-p. [ 2u k R 2 + 0(e\\R\\ 2 2 + 4) + 0(e- 1 \\Rf). 
B k Jn 



Using Lemma loTTl again. this implies 



.4 



A,, 



(6.8) 



— E -E k =- 

B B k 



3E k ' 



±E k f a 2u k R A 



Bi 



TT / 2u k R* 
ok Jn 



0{e\\R\\ 2 2 + + P p )+0{e- l \\R\\ 



The term — ^ J n 2u k R 2 is harmless since it will be handled by the spectrum estimate 
Lemma l6~3l Let us analyze the fraction which is the worst term. For M > 1 to be selected 
later, define || i?|| r , ||-R||t-, the L 2 norms of i? in the tube T := {(x,t) : \d(x,t)\ < Me}, 
the complementary of the tube respectively: 



II?" 2 



T 



f R 2 (x,t)dx, \\R\\ 2 T c := / 

J(\d(x,t)\<Me} Jf\d 



' {\d(x,t)\<Me} 

Observe that the 0(e) size of the tube allows to write 



{\d(x,t)\>Me} 



R 2 (x,t)dx. 



u k R 



T 



< 



T 



1/2 



R- 



T 



1/2 



< Cy/i\\R\\ T . 



Hence, using Lemma loTTl cutting J Q = J r + J 7 - c , we get 

<Ce\E k \(Vi\\R\\T+\\R\\T°) 
<Ce^e\E k \ \\R\\t + Ce 2 / 5 E k 2 + Ce?/ 6 \\R\$ B . 



jT Ek 



As a result 



iEk 2 + j±E k J Q 2u k R > ( 3 _ Ce 2 / 5 )E k 2 - Cesfe\E k \ \\R\\ T 



(6.9) 



Bk 



- Ce 3/5 \\R\\ 2 TC 



B k 



for small e > 0. Now, observe that 
(6.10) E k 2 - Ce^\E k \\\R\\ T > 



if \E k \ > Ce^\\R\\ r 

-C 2 e 3 \\R\\ 2 r if \E k \<Ce^\\R\\T- 



Remark 6.2. The above inequality is the crucial one. One can interpret it as follows. 
Following (8, Proposition 2], we understand that E k behaves like the integral on the hy- 
persurface T t : 

£ / R(x, t) da. 

Jd(x,t)=0 

If \E k \ = \ J n (l - u k 2 )R\ is large w.r.t. 0{e^\\R\\ r ) then E k 2 - Ce^~e\E k \\\R\\ T > 0, 
which has the good sign to control the L 2 norm of R. In other words, if the error "intends" 
at concentrating on the hypersurface, the situation is quite favorable. On the other hand, 
if \E k \ = |/ n (l — Uk 2 )R\ is small w.r.t. C(e^/e||i?||7-) then we get the negative control 
— 0(e 2 1| j-) (after dividing by Bk) which is enough for the Gronwall's argument to work. 
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(6.11) 



Putting together $6J\ . flgjQ , ( l6T9l l, d67TOb and B k = 2^/2\U\e + 0(e 2 ), we arrive at 
A>-^ / 2u fe i? 2 -Ce 3 / 5 ||i?||^-C7 £ 2 ||i?!|2- 



0(e fe+ i)||i?||+0( e || J R||^) + 0(e- 1 || J R|| 3 ). 



6.2. Proof of Theorem l2.3l Equipped with the accurate estimate ( 16.1 U . we can now con- 
clude the proof of the error estimate by following the lines of lfl4l . Combining ( 16.4b with 
(l67TTb and using the interpolation inequality \\R\\l\ P p < C||-R|| p l|Vii|| 2 , \\R\\t < H^H and 
||-R|| ~ 0(e 2 ) (thanks to the definition of t e ), we discover 



i«ii>ii<-/w + ? 



+ ^2d||i?|n|Vi?|| 2 + ^Ce 3 ^\\R\\l e 
+ C\\R\\ 2 + 0(e k -i)\\R\\. 
Since e 2 (- J Q |Vi?| 2 + £ /„(/'(«*) + ^2 Ufe )# 2 ) < -f 2 || Vi?|| 2 + C*||i?|| 2 , we get 



||i?||^||i?|| <(1 - e 2 ) f - |Vi?| 2 + - ^ (f K) + -^2 Ufe ) i? 
-e 2 !|Vi?i| 2 + 42Ci||i?|| p ||Vi?|| 2 
(6.12) + lc7e 3 / 5 ||i?||?- c + C7||i?j| 2 + 0(e fc -^)!|i?|| 

<(l-e 2 ) (-[ |Vi?| 2 + l/ (f'( Uk ) + ^2u k )R 2 



£ Jn\ B k 
+ \ce 3/5 \\R\\?rc + C\\R\\ 2 + 0(£ fc -^)||i?||, 

in view of the definition of t e in (16.5b . In the above inequality, let us write J n = Jj- + JL C . 
In the complementary of the tube, observe that 

/ ( /'(«*) + #2 Ufc + Ce 3/5 ) i? 2 = / (/'K) + 0(e 3/5 )) i? 2 , 

JT C V / J {\d{x,t\>Me} K 7 

is nonpositive if M > is large enough; this follows from the form of the constructed 
u k in Section|5] — roughly speaking we have u k (x,t) = 6q ^ d ( x -*)+ c '( E V j _|_ £)(e 2 ) — 
o (±oo) = ±1 and /'(±1) < 0. As a result we collect 

\R\\j t \\R\\ <(1 - e 2 ) (- ^ |Vi?| 2 + 1 J (/'(u fc ) + ^2u k ^ R 2 



tor —A — £ 2 {f'{u k ) + 4p-2uk) arises when studying the local unbalanced Allen-Cahn 



(6.13) 

+ c7|| J R|| 2 +0(£ fc -^)||i?||. 
In some sense, the problem now reduces to a local estimate since the linearized opera- 

B k 

equation 

d t u e = Au E + 1 (f(u e ) - ^(1 - Ue 2 ) 

whose singular limit is "mean curvature plus a forcing term" (see, among others, J2J). 
To conclude we need a spectrum estimate of the unbalanced linearized operator around the 
approximate solutions u k , namely — A — s~ 2 (f (u k ) + j^2u k ). This directly follows from 
the result of Ifl3l for the balanced case. For related results on the spectrum of linearized 
operators for the Allen-Cahn equation or the Cahn-Hilliard equation, we also refer to Q, 
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Lemma 6.3 (Spectrum of the unbalanced linearized operator around u k lfl3l ). There is 
C* > such that 

- / 1 Vii i 2 + h J (/'(«*) + fi2 ^ c * J r2 i 

for allO <t<T, allO <s<l, all R £ H 1 ^). 
Proof. Observe that 

r 0o (^) +O ( £ 2) if|d (a;>t )|<^i 

u k (x,t) = < V 5 / 

t±l + C(e fe+1 ) if |d(x,t)| > Ve- 

Lemma RTTI yields ^ = e^- + 0(e 2 ) so that we can write f{u k ) + ^2u fc = f'(uk), 
for some 717; such that 

O (^) -e^0i(*J§&) +0(e 2 ) if\d(x,t)\<^ 
±l + 0(e) if\d(x,t)\>y/e, 

where 0! = 1. In particular J R 6»i(rV) 2 /"(#o) = JJPa') 2 f" ( d o) = (odd function) so 
that Ufc has the correct shape for ifPJl to apply: see J4] shape (3.8) and proof of Theorem 
5.1], lfT4l shape (16)] or l22l Section 4] for very related arguments. Details are omitted. □ 

Combining the above lemma and (16.131 l. we end up with 

j t \\R\\<C\\R\\ + Ce k -i. 
The Gronwall's lemma then implies that, for all < t < t £ , 

\\R(;t)\\<(\\R(-M\+e k -he CU =0(e k -h 
in view of ( 16. 11 1. Since k — | > |, this shows that t e = T and that the estimate 0(s 4 / p ) is 
actually improved to 0(e k ~ 5 ). This completes the proof of Theorem l2.3l □ 

Acknowledgements. The authors are grateful to Remi Carles and Giorgio Fusco for valu- 
able discussions on this problem. 
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